Proton-neutron pairing energies in N = Z nuclei at finite temperature 
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Thermal behavior of isoscalar (t=0) and isovector (r=l) proton-neutron (pn) pairing energies 
at finite temperature are investigated by the shell model calculations. These pn pairing energies 
can be estimated by double differences of "thermal" energies which are extended from the double 
differences of binding energies as the indicators of pn pairing energies at zero temperature. We found 
that the delicate balance between isoscalar and isovector pn pairing energies at zero temperature 
disappears at finite temperature. When temperature rises, while the isovector pn pairing energy 
decreases, the isoscalar pn pairing energy rather increases. We discuss also the symmetry energy at 
finite temperature. 
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The proton-neutron (pn) pairing energies have become 
one of hot topics in the study of the nuclear structure for 
proton-rich nuclei. In particular, interests are increas- 
ing in studying isovector (r=l) and isoscalar (r=0) pn 
pairing energies in medium mass N = Z nuclei produced 
at the radioactive nuclear beam facilities. The study of 
pn pairing energies is also important in the astrophysical 
context. These nuclei lie along the explosive rp-process 
nucleosynthesis path and the nuclear properties such as 
masses, halfhves, and isomers have a strong influence on 
modeling the rp-process and identifying possible nucle- 
osynthesis sites. Odd-odd N = Z nuclei are an ideal 
experimental laboratory for the study of pn pairing en- 
ergies. It is well known that the lowest r = and r = 1 
states compete for the ground state changing the sign of 
the energy difference E T —\ — E T= o in odd-odd N = Z 
nuclei, while all even-even N = Z nuclei have the r = 
ground states. Several authors 0, 0, IE El IE IE already 
pointed out that this degeneracy in odd-odd N = Z 
nuclei reflects the delicate balance between the symme- 
try energy and the like-nucleon neutron-neutron (nn) (or 
proton-proton (pp)) pairing energy. On the other hand, 
it has recently been shown that this degeneracy is at- 
tributed to competition between the isoscalar and isovec- 
tor pairing energies [E El HU . 

It has been recently reported [TlllT^ | that the canonical 
heat capacities extracted from observed level densities in 
162 Dy, 166 Er and 172 Yb display the S shape with a peak 
around T « 0.5 MeV, which is interpreted as the break- 
ing of like-nucleon J — pairs because the BCS critical 
temperature corresponds to T c ks 0.57A„(T = 0) w 0.5 
MeV, where the like-nucleon pairing gap A„(T = 0) is 
calculated at zero temperature by the BCS theory. Thus 
it seems that the S shape is a signature of pairing tran- 
sition at the critical temperature. For the finite Fermi 
system like a nucleus, however, since the nuclear radius 
is much smaller than the coherence length of the Cooper- 
pair, statistical fluctuations beyond the mean field in 
the BCS theory become large. The fluctuations smooth 
out the sharp phase transition, and then the like-nucleon 
pairing gap A„ does not quickly become zero at the BCS 



critical temperature but decreases with increasing tem- 
perature. There are many approaches to treat the fluc- 
tuations beyond the mean field. The shell model calcula- 
tion can take into account the large fluctuations beyond 
the mean field. Recently the shell model Monte Carlo 
(SMMC) calculation HEI3 using the fp + g g/2 shell has 
been performed in the even- and odd-mass Fe isotopes. 

We recently proposed "thermal" odd-even mass 




FIG. 1: r = and r = 1 pn pairing gaps estimated from dou- 
ble differences of binding energies for odd-odd N = Z nuclei: 
(a) experimental ones; (b) those of shell model calculations. 
The solid and dotted curves show 122.25(1 - 1.67A~ 1/3 )/A 
and 5.18A -1 / 3 , respectively. 
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difference to estimate the like-nucleon pairing energy at 
finite temperature, and showed in the spherical shell 
model calculations that the drastic suppression of like- 
nucleon pairing energy due to finite temperature brings 
about the S shape in the heat capacity arround the tem- 
perature T c w 0.57A„(T = 0) MeV. In this rapid com- 
munication, we study the pn pairing energies at finite 
temperature in odd-odd N = Z nuclei. Does pairing 
transition due to the breaking of pn pairs take place 
when temperature increasing? It is now interesting to 
investigate thermal behavior of the pn pairing energies 
in N = Z nuclei. 

We start from the double difference of binding energies 
[3 13 El defined as 



& T pn (Z, N) = -[B(Z, Nf - B{Z, N-l) 
-B{Z — 1,N) + B(Z — 1,N— I)], 



(1) 



where B(Z, N) is the binding energy. The indicator A£„ 
gives the r = 1 pn pairing gap in N = Z nuclei. The 
AJ^ can be regarded as the r = pn pairing gap as 
well. Figure 1 (a) shows the r = and r = 1 pn pair- 
ing gaps estimated from the double differences of exper- 
imental binding energies (I) in odd-odd N — Z nuclei 
with A — 18 — 58. The r = pn energy is somewhat 
larger than the r = 1 pn energy in the sd shell nuclei and 
vice versa in the pf shell nuclei. Over a wide range of 
odd-odd N — Z nuclei, however, basically Fig. 1 shows 
almost the same magnitude of the t = and r = I pn 
pairing gaps. We carried out shell model calculations 
using isospin-invariant interactions such as the unified 
sd (USD) interaction E3 f° r odd-odd N = Z nuclei in 
sd shell and the GPFX1 interaction [H for 42 Sc, 46 V, 
and 50 Mn in fp shell. On the mean-field level the ra- 
tio between the strengths of pp-, nn-, and pra-pair fields 
is given by the orientation of the pair field. The rela- 
tive strengths of three types of pair fields becomes only 
definite when isospin symmetry is restored. Note that 
the shell model calculations with isospin invariance show 
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Ap" 1 in odd-odd N = Z nuclei. 



In Fig. 1(b), we can see that the shell model results 
reproduce well the experimental pn pairing gaps, and de- 
scribe the characteristic behavior in Fig. 1(a). The pn 
pairing gaps are closely realted to the energy difference 
B(Z, N) T=1 - B(Z, iV) r=0 between the lowest r = and 
r = 1 states in odd-odd N — Z nuclei, because the en- 
ergy difference satisfies the following identity 0, 



B(Z, N) 7 



B(Z,N) 



r=0 



T=0 



A 



T = l\ 



(2) 



Odd-odd N = Z nuclei with A < 40 have the ground 
states with r = 0, J > except for 34 C1, while the 
ground states of odd-odd N = Z nuclei with 40 < A < 74 
are r = 1 and J = except for 58 Cu. Several authors 
0i S IE discussed that this degeneracy is at- 
tributed to the delicate balance between the symmetry 
energy cl(A)t(t + 1) /A and pairing gap A and that the 
energy difference SB = B(Z, N) T=1 - B{Z, N) T =° is ex- 



pressed as SB = 2(a(A)/A — A). However, if we em- 
ploy the symmetry energy coefficient a(A) = 134.4(1 — 
1.52 A' 1 / 3 ) and pairinggap A = 5A8A- 1 / 3 of Dufio and 
Zuker mass formula |2l|. the energy difference SB be- 
comes larger than the experimental value. As suggested 
in our previous paper, the isoscalar pairing gap A^ 

is approximately written as 122.25(1 - 1.67A~ 1/3 )/A 
and the isovector one AJ^ 1 is equal to the like-nucleon 
nn pairing gap A„ « 5.18A -1 / 3 . These two curves 
are shown in Fig. 1(b) for comparison. Since SB = 
2(Ap J " 1 — Ap=°), the degeneracy between the lowest t = 
and t = 1 states in odd-odd N — Z nuclei comes from 
the delicate balance between the isoscalar and isovector 
pn pairing energies. 

Let us next describe the pn pairing gaps at finite tem- 
perature. We introduce the canonical partition function 
defined by 



Z(T) = Tr(e-»/ T )=J2< 



-Ei/T 



(3) 



i=0 



where Ei is the energy of the ith eigenstate with degen- 
eracies based on symmetries for the Hamiltonian H of a 
system. All the eigenvalues Ei are obtained by solving 
the eigenvalue equations H^i — E^i. Then, the parti- 
tion function in the canonical ensemble is calculated from 
Eq. J3J, and any thermodynamical quantities 0(T) can 
be evaluated from 



0(T) = (O) = Tr(Oe- ff / T )/Z(T), 



(4) 



where (O) stands for the average value of operator O 
over the range of eigenstates. For instance, the thermal 
energy is expressed as 



E(Z,N,T) = (H) = Y d E i e~ E ^ T /Z{T). 



i=0 



The heat capacity is then given by 

dE(Z, N, T) 



C{Z,N,T) 



dT 



(5) 



(6) 



We now introduce the following double difference of 
"thermal" enegies E(Z, N,T) analogous to Eq. (1) as an 
indicator of pn pairing energies, 

A; n (Z, N, T) = l -[E(Z, N, T) T - E(Z, N - 1, T) 

-E(Z-1,N,T) + E(Z -1,N -1,T)]. (7) 

The double differences of binding energies at zero tem- 
perature in Eq. (1) are known theoretically and experi- 
mentally as important quantities in evaluation of the pn 
pairing energies in a nucleus [id E1E3- The double dif- 
ferences of thermal energies in Eq. (7) are also indicators 
of the pn pairing energies and can be regarded as the pn 
pairing gaps at finite temperature. 

Let us evaluate the double difference of thermal ener- 
gies (7) for N — Z sd shell nuclei. We make numerical 



3 



calculations by way of two steps. First, we carry out the 
exact shell model calculations in the sd shell using the 
USD interaction {B^, and calculate the correlated ther- 
mal energy E v>tr from Eq. JSJ. Secondly, we extend the 
model space to a larger one (sd + fp + s 1 / 2 d 5 / 2 ) in or- 
der to display the double difference of thermal energies 
in a broader range of temperature using an independent- 
particle approximation |14|. The single-particle energies 
of the extended space are obtained by diagonalizing the 
Woods-Saxon potential with the spin-orbit interaction, 
where the harmonic-oscillator (H.O.) cigenfunctions are 
used. The Woods-Saxon parameters are chosen so as 
to reproduce the single-particle energies estimated from 
17 0, because it is necessary to reasonably extrapolate 
the single-particle energies of the sd shell to those of the 
larger space. In this way, we combine the correlated ther- 
mal energy E Vttr in the truncated space with the ther- 
mal energy E sp calculated using the independent-particle 
approximation in the larger space. The thermal energy 
which takes account of the interaction effects in the sd 
shell is estimated as follows UM: 



E 



E v tr + E ap — E s 



(8) 



where E sp ^ r is the thermal energy of the sd shell within 
the independent-particle approximation. We now obtain 
the double difference of thermal energies A pn by substi- 
tuting E of Eq. © for E(Z, N, T) in Eq. (7). 

Figure 2 shows the calculated thermal pn pairing gaps 
for odd-odd N = Z nuclei, 22 Na, 26 Al, 30 P, and 34 C1 at 
temperature T — 2.0 MeV. The r =1 and r =0 pn pairing 
gaps are largely separated at T — 2.0 MeV. Comparing 
Fig. 2 with Fig. 1(b), we notice that the r =1 pn pairing 
gap decreases but the r =0 pn pairing gap keeps the 
magnitude from zero temperature to high temperature. 

Figure 3 shows the variation of the thermal pn pairing 
gaps depending on temperature T for 22 Na, 26 Al, 30 P, 
and 34 CI. In all graphs, we can see increase of the r =0 
pn pairing gap and decrease of the t =1 pn pairing gap. 
As mentioned above, at zero temperature the t =0 and 
t =1 pn pairing gaps are almost the same, and the lowest 
t =0 and t =1 states are degenerate. As increasing 
temperature, the r =1 pn pairing gap decreases and r =0 
one rather increases. Thus, we know that the r = 
pairing energy becomes dominant at high temperature. 

It would be valuable to discuss the symmetry energy 
~ 4a sym (T)T(r + 1)/A at finite temperature because it is 
closely related to the r = pairing energy. In our pre- 
vious paper [l^, we suggested that the dominant part 
of the symmetry energy comes from the r = pairing 
energy part in the shell-model interaction energy. For 
the application of the symmetry ener gy i n core-collapse 
supernova simulations, Donati et al. [23j pointed out a 
possibility that the symmetry energy coefficient a sym at 
the finite temperature has been estimated to be some- 
what larger than that of stable nuclei at zero tempera- 
ture. The increase (~ 3%) of the symmetry energy be- 
tween T — 0.0 and T = 1.0 MeV after implementing the 
correction in the SMMC calculations is smaller than that 




FIG. 2: Calculated thermal pn pairing gaps for odd-odd N = 
Z nuclei at temperature T = 2.0 MeV. The solid circles denote 
the t—0 pn pairing gap, and the open circles the r=l pn 
pairing gap. 



(~ 8%) of the quasiparticle random phase approximation 
(QRPA) with temperature |24| . We now calculate the 
temperature dependence of the symmetry energy using 
the shell model calculations. We estimate the symmetry 
energy coefficient from the thermal energy E(Z, N,T) T 
with isospin r and temperature T obtained in the shell 
model calculations as follows: 



E(Z, V,T) r - E(Z, N, T) T 
t(t + 1)-t'(t' + 1) 



A, (9) 



where r and r' are different isospins for isobaric nuclei 
with same mass number A. At zero temperature, the 
calculated symmetry energy coefficient a sym (T = 0) ~ 16 
MeV for A=2A is in good agreement with the value deter- 
mined from experimental masses and with the empirical 
value of Duflo and Zuker mass formula. 

Figure 4 shows the symmetry energy coefficient a sym 
as a function of the temperature for even-even N f» Z 
nuclei with mass number A =20, 24, and 28, where sev- 
eral isobaric pairs of N ps Z nuclei such as ( 20 Ne, 20 O), 
( 24 Mg, 24 Ne), and ( 28 Si, 28 Mg) are chosen. This figure 
shows that the symmetry energy coefficients increase 
with increasing temperature in these three cases. More- 
over, we can see that the symmetry energy coefficient 
depends on the mass A which is empirically fitted by 
adding the surface contribution with the A~ 4 / 3 depen- 
dence at zero temperature. This mass dependence ap- 
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FIG. 3: Calculated thermal pn pairing gaps for odd-odd 
N — Z nuclei as a function of temperature. The solid line 
denotes the r=0 pn pairing gap, and the dotted line the t=1 
pn pairing gap. 



pears in the r = pn pairing gap estimated from the 
double difference of binding energies, in Fig. 1(b). Fig- 
ure ^] also suggests that the mass dependence changes as 
temperature increases. To see the temperature depen- 
dence of the symmetry energy coefficient, we define the 
relative change of the symmetry energy coefficient with 
respect to temperature as 



Sa sym (T) = 



' syni 



(T) - a sym (T = 0) 



a sym (T = 0) 



(10) 



Averaging the Sa synl (T) at T = 1.0 MeV over various 
pairs of nuclei, we obtain an increase ~ 4%. This is in 
agreement with the SMMC result ~ 3% obtained after 
implementing the correction. We used here the form of 
symmetry energy r(r + 1), which is motivated by the 
charge independence of the nuclear force. But as a phe- 
nomenological parametrization the isospin dependence 
t(t + a) with a ^ 1 is also possible, where the linear 
term in r is so-called Wigner term. Recently, empirical 
fitting to the Wigner ter m g ave a = 1.25 in the vicin- 
ity of the N = Z line 0, H3. However, the symmetry 
energy coefficient is affected little by replacing r(r + 1) 
with t(t + 1.25). Moreover, by definition the relative 



change of the symmetry energy coefficient <5a sym (T) does 
not change by this replacement. 

In conclusion, we investigated the t—0 and r=l pn 
pairing energies at finite temperature using the shell 
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FIG. 4: Symmetry energy coefficient a sym as a function of 



temperature for N 
and 28. 



Z nuclei with mass number A =20, 24, 



model calculations. The pn pairing gaps at finite tem- 
perature were estimated from the double differences of 
thermal energies defined by Eq. (7), which is analogous 
to the double differences of binding energies as indica- 
tors of the pn pairing energies at zero temperature. It 
was shown that as temperature increases the isoscalar 
pn pairing energy rather increases, while the isovector 
pn pairing energy decreases. Almost the same pn pair- 
ing gaps of stable N = Z nuclei at zero temperature are 
separated with increasing temperature. We also stud- 
ied the temperature dependence of the symmetry energy 
in N rj Z nuclei. The symmetry energy coefficients in- 
crease with increasing temperature. The increase of the 
calculated symmtery energy coefficient between T = 0.0 
and T = 1.0 MeV is in good agreement with that of the 
SMMC calculations. We suggest that the pn pairing en- 
ergies can be estimated using Eqs. (5) and (7) from the 
measured level densities of nuclei. We expect that the 
pn pairing energies play an important role in the astro- 
physics. 
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